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ABSTRACT 

We have analyzed seventeen observations of the low-mass X-ray binary and atoll 
source 4U 1728-34, carried out by the Rossi X-ray Timing Explorer in 1996 and 1997. 
We obtain precise measurements of the frequencies of the two simultaneous kilohertz 
quasi-periodic oscillations (kHz QPOs) in this source. We show that the frequency 
separation between the two QPO, Au, is always significantly smaller than the frequency 
of the nearly-coherent oscillations seen in this source during X-ray bursts, even at the 
lowest inferred mass accretion rate, when Av seems to reach its maximum value. We 
also find that Au decreases significantly, from 349.3 ± 1.7 Hz to 278.7 ±11.6 Hz, as the 
frequency of the lower frequency kHz QPO increases from 615 to 895 Hz. This is the first 
time that variations of the kHz QPO peak separation are measured in a source which 
shows nearly-coherent oscillations during bursts. 

Subject headings: accretion, accretion disks — stars: neutron — stars: individual (4U 
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1. Introduction 

Quasi-periodic oscillations at frequencies near 1 
kilohertz (kHz QPOs) are now known to be a com- 
mon phenomenon among low mass X-ray binaries 
(LMXBs). Since the beginning of the Rossi X-ray 
Timing Explorer (RXTE) mission kHz QPOs have 
been discovered in the persistent flux of nineteen 
LMXBs (see van der Klis 1998| for a review). In 
almost all cases the power density spectra of these 
sources show twin kHz peaks that, as a function of 
time, gradually move up and down in frequency, typ- 
ically over a range of several hundred Hz. 

The dependence of the QPO frequencies upon X- 
ray flux is compl ex. One example is 4U 1608-52 
(Mendez et al. 1999| ): While on time scales of 
hours frequency and X-ray flux are well correlated, 
at epochs separated by days to weeks the QPOs span 
the same range of frequencies even if the average flux 
is different by a factor of 3 or more (see also Ford 
et al. [19971; Zhang et al. |1998|). Ford et al. (|l997j) 



have shown that in 4U 0614+09 the QPO frequency is 
better correlated to the flux of the soft ( "blackbody" ) 
component in its X-ray spectrum than to the total X- 
ray flux, while Kaaret et al. (1998) found, both for 
4U 0614+09 and 4U 1608-52, a good correlation be- 
tween the QPO frequency and the photon index of the 
power- law co mponent in the X-ray spectrum. Mendez 
et al. ( 1999 ) showed that QPO frequency is very well 
correlated to the position of the source on the color- 
color diagram, and they concluded that most likely 
these frequency changes are driven by changes in M 
through the innermost part of the accretion disk, and 
that there is no one-to-one relationship between the 
observed X-ray flux and this mass accretion rate. 

Initially, data from various sources suggested that 
the separation Av between the twin kHz peaks re- 
mained constant even as the peaks moved up and 
down in frequency (e.g., Strohmayer et al. 1996] ). In 
some sources a third, nearly-coherent, oscillation is 
detected during type-I X-ray bursts, at a frequency 
close to Av, or twice that value (see Strohmayer, 
Swank, & Zhang 199S for a review). These two re- 
sults suggested that a beat-frequency mechanism was 
at work (Strohmayer et al. 1996| ; Miller, Lamb, & 
Psaltis 1998j ), with the third peak being close to the 
neutron star spin frequency or twice that. 

It is now clear that Ais, which in the beat-frequency 
interpretation is thought to be equal to the spin fre- 
quency of the neutron star, is not always constant. 



In Sco X-l (va n der Klis et al. [1997D , 4U 1608-52 
(Men dez et al. |l998c[ ), and 4U 1735-44 (Ford et al. 
1998| ) , Ai> slowly decreases as the frequencies of both 
QPOs increase. In these sources no burst oscillations 
have been detected so far. In 4U 1636-53 the fre- 
quency of the burst oscillations (or half this value) 
does not exactly match the frequency difference be- 
tween the kHz QPOs (Mendez, van der Klis, & van 
Paradijs 1998a). In this case, however, no evidence 
was found of any significant change of Av. So far, 
no measurements existed of burst oscillations and a 
varying kHz peak separation in the same source. 

In this paper we present such measurements for the 
first time. We show that in 4U 1728-34, the frequency 
separation between the two simultaneous kHz QPOs 
is always lower than the frequency of the burst oscil- 
lations, even at low kHz QPO frequencies, when the 
peak separation appears to have reached its highest 
value. We also find that Av decreases significantly as 
the frequencies of the kHz QPOs increase. As in 4U 
1608-52, the QPO frequencies in 4U 1728-34 are well 
correlated to the position of the source on the color- 
color diagram, despite them not being correlated to 
the X-ray flux. 

2. Observations 

RXTE observed 4U 1728-34 on ten occasions in 
1996, between February 15 and March 1, for a total 
of 199 ks (for a de script i on of these observations see 
199(1 " 



Strohmayer et al. 



1997: Ford & van der Klis 



1998| ), and seven more times starting UTC 1996 May 



3 13:56:00 for 3.9 ks, 1997 September 23 23:48:00 for 
1.4 ks, September 24 09:13:00 for 22.9 ks, September 
26 12:28:00 for 18.1 ks, September 27 09:17:00 for 18.3 
ks, September 30 04:31:00 for 20.6 ks, and October 1 
06:07:00 for 21.2 ks. 

Here we only use data recorded by the Proportional 
Counter Array (P CA) onboard RXTE (Bradt, Roth- 
schild, & Swank 1993] ). Besides the two standard 
modes that are available to all RXTE/PCA obser- 
vations, data were collected using additional modes 
with high time, and moderate energy resolution, cov- 
ering the nominal 2 — 60 keV energy band of the 
RXTE/PCA. Part of these data (the 1996 observa- 
tions) were analyzed previously (Strohmayer et al. 
1996| , [1997D - Their analysis of the kHz QPOs was 
based on a simple direct FFT method, while ours uses 
the "shift-and-add" method (see below). 

To estimate the source X-ray intensity and colors 
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we measured the count rates every 16 s in 5 energy 
bands, 2.0 - 3.5 - 6.4 - 9.7 - 16.0 keV, and 2.0 - 16.0 
keV, taking into account the gain changes applied to 
the PCA in March and April 1996. In a few of the 
observations one or two of the five detectors of the 
PCA were switched off; we only used the three de- 
tectors which were always on to calculate these count 
rates. We subtracted the background contribution in 
each band using the standard PCA background model 
version 2.1b, and normalized the count rates to 5 de- 
tectors. For the timing analysis we used data from all 
active detectors without energy selection. 

3. Results 

3.1. Measurement of the frequencies of the 
QPOs 

We divided the high-time resolution data into seg- 
ments of 64 s, and calculated a power spectrum 
for each segment preserving the maximum available 
Nyquist frequency (8192 or 65536 Hz, depending on 
the observation). We produced an average power 
spectrum for each observation and searched it for 
QPOs at frequencies above 100 Hz. 

In the rest of this subsection we only present the 
results of the analysis of those 1482 data segments in 
observations where we detected two simultaneous kHz 
QPOs in the power spectrum. Here we concentrate 
on the measurements of the frequency separation Av 
of these QPOs, as other properties of the QPOs have 
already been reported in previous papers (Strohmayer 
et al. 199£ ; Ford & van der Klis 1998; for the observa- 



tions that we report here for the first time, the other 
properties of the QPOs are similar to those described 
in these two papers). 

To measure Av we applied the 'shift-and-add' tech- 
nique advanced by Mendez et al. (1998b). Measur- 
ing Av in an averaged (non-shifted) power spectrum 
can introduce a bias that changes the peak separation 
from its true value because, in general, the two QPOs 
vary in strength in a different manner as the y mov e 
in frequency (see for instance Wijnands et al. 1997 ). 
This is no longer the case if the individual power spec- 
tra are shifted to the frequency of one of the two 
peaks before they are averaged together. (However, 
this method can be applied only if one of the QPOs is 
always strong enough to measure its frequency in each 
of the individual power spectra.) A second advantage 
of this method is that longer stretches of data can be 
used to measure Av, and not just pieces of data where 



the QPO frequencies are more or less constant. This, 
in turn, decreases the errors in the measurements of 
the frequency separation (see below). 

We measured the central frequency of the lower 
frequency QPO in each segment of 64 s by fitting a 
Lorentzian to it (this QPO was well detected in each 
segment). We then grouped the data in 8 sets of 30, 
58, 66, 193, 371, 262, 244, and 258 power spectra, such 
that the frequency of the QPO did not vary by more 
than ~ 50 Hz within each set (the first 2 sets contain 
data from 1996 only, whereas the last 2 sets contain 
data from 1997 only; segments 3 to 6 include data 
from both epochs). We then shifted the frequency 
scale of each spectrum to a frame of reference where 
the position of this peak was constant in time, and 
we finally averaged the aligned spectra in each set. 

We fitted the resulting 8 average power spectra in 
the range 256 — 1500 Hz using a function consisting 
of a constant, representing the Poisson noise, and two 
Lorentzians, representing the QPOs. The fits were 
good, with reduced x 2 < 1-1 f° r 1585 degrees of free- 
dom, and the significance of each peak was always 
> 4cr. In Figure [j] we plot the frequency difference, 
Av, between the higher frequency and the lower fre- 
quency QPOs (hereafter the upper and lower QPO, 
respectively) as a function of the centroid frequency of 
the lower QPO, v\ ow . The solid line indicates the fre- 
quency of the oscillations detected during the X-ray 
bursts, t'burst = 363.95 Hz (Strohmayer et al. 1996, 
1998| ) . These results are entirely consistent with those 
of Strohmayer et al. ( |1996| ; Table 1), but our error 



bars, thanks to the use of the "shift-and-add" tech- 
nique, are much smaller than theirs. 

Because of these small error bars, it is now ap- 
parent that Av is always smaller than t'tmrst- The 
figure shows that Av is consistent with being con- 
stant at v\ OVJ < 800 Hz, but decreases as v\ ow increases 
above 800 Hz. The data point near v\ ow — 900 Hz 
is significantly (6a) below the average of the other 
points. If, conservatively, we discard both points at 
^low > 800 Hz the average frequency separation be- 
comes Av = 349.3 ± 1.7 Hz (indicated by the dashed 
line in Fig. |l|), still different from 363.95 Hz at an 
8.6cr level. So, even in the range where Av appears 
to have reached its saturation value, it is significantly 
below the frequency of the burst oscillations. 
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3.2. QPO frequency as a function of inferred 

M 

In Figure || we show a color-color diagram of 4U 
1728-34, which is shown here for the first time to 
be typical of t hose o f other atoll sources (Hasinger 
& van der Klis 1989 ). In general, the count rate is 
observed to increase as indicated by the arrow on the 
diagram, from the island state (at the upper right part 
of the diagram) to the banana branch, but the relation 
between count rate and either of the two colors is 
much less clean than that between colors. Each point 
in this diagram represents 128 s of data. Dots and 
open circles indicate the segments with and without 
detectable kHz QPOs, respectively. 

In Figure |^ we show the frequency of the lower 
QPO, J4owi as a function of count rate for 4U 1728- 
34. Each point in this figure represents 128 s of data. 
The branch in the lower right part of the diagram 
corresponds to the data presented by Strohmayer et 
al. ( |1996| ) for the lower kHz QPO. This plot c an be 
compared to those in Aql X-l (Zhang et al. 1998) 
and 4U 1608-52 (Mendez et al. |1999|). As in those 



sources, for 4U 1728-34 it is also true that on time 
scales greater than ~ 1 day there is no unique relation 
between count rate and v\ ow - 

As in the case of 4U 1608-52 (Mendez et al 



1999) 



QPO frequencies correlate much better to the posi- 
tion of the source on the color-color diagram than to 
count rate. In Figure || we show v\ ow as a function of 
the position on the color-color diagram as measured 
by the parameter S a (curve length in the color-color 
plane measured in t he dir ection of the arrow in Fig. 
U see Mendez et al. 1999 for the precise definition of 
Sa), for the same intervals shown in Figure | (black 
circles). In this case S a = 1 at colors (0.48,-0.23), 
and S a = 2 at (0.44, —0.34). The black squares repre- 
sent ^ U pp) the frequency of the upper kHz QPO as a 
function of S a (we only include here segments where 
we could measure ^ upp directly). The grey triangles 
are measurements from segments in which we only 
detect one of the kHz QPOs; however, from the loca- 
tion of each point in this diagram we can determine 
whether it is the lower or the upper QPO. While the 
frequency vs. count rate relation (Fig. ||) is complex, 
there seems to be a neat relation between each QPO 
frequency and S a (Fig. ^). 



4. Discussion 

We have for the first time measured the variation 
of the kHz QPO peak separation in a source which 
shows oscillations during X-ray bursts. The result of 
this measurement is unexpected: even at its satura- 
tion level at low inferred M, Av is still significantly 
below the frequency of the burst oscillations. Before 
we discuss this further, we consider what our data 
imply about the relation between the kHz QPO fre- 
quencies and M. 

Despite the complex relation between frequency 
and count rate (Fig. ||), in 4U 1728-34 there is a 
unique relation between frequency and position in the 
X-ray color-color diagram (Fig. ^), similar to what 
is observed in 4U 1608-52 (see Fig. 2, 3, and 4 in 
Mendez et al. |l999t ). As in 4U 1608-52, in 4U 1728- 
34 the frequencies of the QPOs increase with S a , as 
the source moves from the island to the banana. The 
presence of the QPOs also correlates with the position 
in the color-color diagram. This may be a general 
feature. In atoll sources it has been proposed that 
M increases with S a from the island to the banana 
(Hasinger & van der Klis 1989), even though X-ray 
flux is not monotonically related to S a (van der Klis 
et al. 119901; van der Klis 11994!; Prins & van der Klis 



1997). Our proposed interpretation is that the ob- 
served changes in the frequency of the kHz QPOs in 
LMXBs are driven by changes in M through the in- 
nermost part of the accretion disk, and that while 
the observed X-ray flux is not a good indicator of this 
mass accretion rate, the X-ray colors track it much 
better. Possible reasons for this include anisotropies 
in the emission, redistribution of some of the radia- 
tion over unobserved energy ranges, and matter that 
flows onto the neutron star in a more radial inflow, or 
away from it in a jet. 

We have found that in 4U 1728-34 the frequency 
separation between the two simultaneous kHz QPOs, 
Az/, is not consistent with being equal to the fre- 
quency of the burst oscillations. This is the sec- 
ond source, after 4U 1636-53, for which Av does 
not match i^burst- Furthermore, in 4U 1728-34 Av is 
not constant, but it decreases as the frequency of the 
lower kHz QPO increases. This is the first source in 
which the three oscillations are seen simultaneously, 
and in which the separation between the frequencies 
of the kHz QPOs is observed to change. 

It has been proposed that the upper kHz QPO 
at f upp represents the Keplerian frequency of the ac- 
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creting material in orbit around the neutron star at 
some preferred radius (van der Klis et al. 1996| ), while 
the lower QPO, at v\ ow , is produced by the beating 
of JAipp with the spin frequency of the neutron star, 
^spin, which has been identified with Z4> ur st (or half 



that; Strohmayer et al. 1996; Miller et al. 1998) 



If this interpretation is correct, Av equals v sp [ n and 
should remain constant. The neutron star spin can 
not change by ~ 20 — 30 % on timescales of days to 
months. 

There are now four sources for which Av has been 
observed to vary: Sco X-l (van der Klis et al 



1997) 



4U 1608 -52 ( Mendez et al. pMcj , 4U 1735-44 (Ford 
et al. p9"8| ), and 4U 1728-34. It has been ar- 
gued that these variations can be attributed to the 
near-E ddington mass accretion rate: White & Zhang 
( 1997 ) propose a 35% expansion of the neutron star 
photosphere with conservation of angular momentum, 
and F. K. Lamb (1996, private communication) sug- 
gests that the height of the inner disk increases, and 
the different values of Av = v upp — i^ ow reflect differ- 
ent values of v upp at different heights in the disk. 

Contrary to what we observe in 4U 1728-34, these 
explanations imply that Av should approach v sp - ln 
when, at low M, the kHz QPOs move to lower fre- 
quencies^ Our results show that in 4U 1728-34, Au 



is significantly below ^burst at the lowest inferred 



M values. This seems to rule out the simple beat- 
frequency interpretation of the kHz QPOs in LMXBs, 
and some of the modifications introduced to explain 
the results of a variable Av. 

Recent ideas allowing the observed beat-frequency 
to deviate slightly from the true one (F. K. Lamb 
1998, private communication), as well as proposals 
discarding one or more elements of the basic beat- 
frequency picture, but still predicting definite rela- 
tions between the observed frequencies (e.g., Stella 



& Vietri |199j; Titarchuk, Lapidus, & Muslimov |199J) 
therefore merit careful consideration. 
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Fig. 1. — The frequency separation between the upper and the lower QPO peak as a function of the frequency of 
the lower QPO for 4U 1728-34 (circles). The solid line indicates the frequency of the burst oscillations at 363.95 
Hz (Strohmayer et a. 1996, 1998). The dashed line indicates the average value of Az/ (349.3 Hz) for the first 6 
points (see text). 
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Fig. 2. — Color-color diagram of 4U 1728-34. The soft and hard colors are defined as (the logarithm of) the ratio 
of count rates in the bands 3.5 — 6.4 keV and 2.0 — 3.5 keV, and 9.7 — 16.0 keV and 6.4 — 9.7 keV, respectively. 
The contribution of the background has been subtracted, but no dead-time correction was applied to the data (in 
this case the dead-time effects on the colors are less than 1 %). Each point represents 128 s of data. Typical error 
bars are shown. Dots and open circles indicate segments with and without detectable kHz QPOs, respectively. The 
arrow indicates the sense in which M is inferred to increase. 
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Fig. 3. — Relation between the frequency of the lower kHz QPO and the 2 — 16 keV count rate. The count rates 
have been corrected for background, and normalized to 5 detectors. Each point represents a 128-s segment. We 
only includes data where both kHz QPOs were detected simultaneously, so that we can unambiguously identify the 
lower kHz peak. 
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Fig. 4. — Diagram of the frequencies of the upper and lower kHz QPOs vs. the position of the source on the 
color-color diagram, as measured by S a (see Fig. |], and text). Black circles and squares represent the lower and 
the upper kHz peak, respectively. Grey triangles represent segments where we only detected one of the kHz QPOs 
in the power spectrum. 
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